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ABSTRACT 

The results of a multisite photometric campaign of 9^ Tau are reported. This bi- 
nary system consists of evolved and main-sequence A stars inside the instability strip. 
The 12th Delta Scuti Network campaign included 152 hours of high-precision photom- 
etry obtained at four observatories. This leads to the derivation of 11 frequencies of 
pulsation in the 10.8 to 14.6 cd~^ range. These frequencies confirm the results from 
previous Earth-based (1982-1986) as well as satellite (2000) photometry, although 
amplitude variability on a time scale of several years is present. 

We show that at least two high frequencies (26.18 and 26.73 cd~^) are also present 
in O'^ Tau. Arguments are given that these high frequencies originate in the main- 
sequence companion and are not combination frequencies, f^-l-fj, from the primary. 

Models for both the primary and the secondary components were checked for the 
instability against radial and nonradial oscillations. All hot models of the primary with 
Teff > 8000 K are stable in the observed frequency range. The best fit between the 
theoretical and observed frequency ranges is achieved for models with Teg w 7800 K 
(or slightly higher), in agreement with photometric calibrations. The instability range 
spans two or three radial orders in the range p4 to pq for radial modes. Post-main- 
sequence models (with or without overshooting) are preferable for the primary, but 
main-sequence models with overshooting cannot be excluded. For the less luminous 
secondary component the instability range is wider and spans 5 to 7 radial orders 
from p2 to ps- The observed frequencies lie around radial modes p^ to pg- The main 
uncertainties of these results are caused by a simple treatment of the convective flux 
in the hydrogen ionization zone. 

Key words: 6 Scuti - Stars: oscillations - Stars: individual: 9^ Tau - Stars: individual: 
78 Tau - Techniques: photometric 



1 INTRODUCTION 

A number of lengthy observing campaigns covering individ- 
ual 5 Scuti variables have shown that the majority of these 
pulsating variables on and near the main sequence pulsate 
with a large number of simultaneously excited nonradial p 
modes. Furthermore, long-term variability of the pulsation 
amplitudes of these nonradial modes with a time scale of 
years has been discovered for many, but not all of these non- 
radial pulsators. The Delta Scuti Network, hereafter called 
DSN, specializes in multisite observations of these stars. The 
network is a collaboration of astronomers located at obser- 
vatories spaced around the globe in order to avoid regular. 



daily observing gaps. So far, 21 campaigns have been car- 
ried out. The most recent campaign covered BI CMi, for 
which 29 frequencies of pulsation derived from 1024 hours 
(177 nights) of photometry were discovered (Breger et al. 
2002). 

The variability of the star 9^ Tau was discovered by 
Koran (1977, 1979) and confirmed by Duerbeck (1978). In 
order to study the multiple frequencies, the Delta Scuti Net- 
work undertook two campaigns (Breger et al. 1987, 1989: the 
latter is referred to as Paper I). The data could be supple- 
mented by additional measurements by Kovacs & Paparo 
(1989). Five frequencies of pulsation were detected and in- 
terpreted to be due to nonradial pulsation. The size and dis- 
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tribution of the frequencies suggested p modes with values of 
£ = to 2. A mode with a higher £ value was also detected 
spectroscopically by Kennelly & Walker (1996). Since the 
original analysis of the photometric data, new theoretical as 
well as experimentally determined (Brogor ct al. 1993) sta- 
tistical criteria for the acceptance or rejection of additional 
modes extracted from multisite photometric data have be- 
come available. A re-examination of the previous data ob- 
tained in the years 1982 to 1986 (Breger 2002, hereafter 
referred to as Paper II) led to the discovery of ten frequen- 
cies in the 1986 DSN data. Nine of these frequencies are also 
present in the list of twelve frequencies determined during a 
remarkable study of 0^ Tau from space (Porctti ct al. 2002, 
hereafter referred to as Paper III): during 2000 August, 9^ 
Tau was monitored extensively with the star tracker camera 
on the Wide-Field Infrared Explorer satellite (WIRE). 

6^ Tau is a 140.728d binary system with known orbital 
elements (Ebbighausen 1959, Torres, Stefanik & Latham 
1997). The two components have similar temperatures. The 
primary component is evolved (A7III), while the secondary 
is fainter by 1.10 mag (Peterson, Stefanik & Latham 1993) 
and still on the main sequence. Both stars are situated in- 
side the instability strip. It was shown in Paper I that the 
dominant pulsation modes originate in the primary compo- 
nent because the predicted orbital light-time effects for the 
primary match the observed shifts in the light curves of up 
to several minutes. This confirms an expectation from the 
values of the frequencies (10 - 15 cd^^), which are incom- 
patible with those of main-sequence 5 Scuti stars. 
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2 NEW GROUND-BASED MEASUREMENTS 

In order to eliminate the serious alieising caused by regular 
observing gaps, a multisite campaign was organized utilizing 
the Delta Scuti Network (DSN campaign 12). During 1994 
November and December, 6^ Tau was measured photometri- 
cally with the Three-Star- Technique (Breger 1993) at a num- 
ber of observatories spaced around the globe. Weather con- 
ditions, unfortunately, wore unfavorable so that data from a 
number of Asian Observatories were not accurate enough to 
be included in the present analysis. 152 hours of photometry 
at four other observatories could be used. These observato- 
ries were: 

(i) Lowell: The photoelectric photometer attached to 
the 0.8m reflector at Lowell Observatory, Flagstaff, USA, 
was used with the Stromgren v and y filters together with a 
neutral density filter. The observer was W. Zima. 

(ii) McDonald Observatory: The 0.9m telescope was 
used with the y filter. The observer was G. Handler. 

(iii) Sierra Nevada Observatory: The 0.9m telescope was 
used together with v and y filters. The observers were R. 
Garrido and F. Boichbuchncr. 

(iv) The 0.6m reflector at the XingLong station of Bei- 
jing Astronomical Observatory, China, was used for 9 nights 
with the V filter. The observers were Li Zhiping, Zhou Aiy- 
ing, and Yang Dawei. More details can be found in Li et 
al. (1997a), whore the Chinese data have already been pub- 
lished, and Li et al. (1997b). We have adopted the Chinese 
data as published except for the slight editing of a few (ob- 
viously) deviant points. 



Figure 1. Light curves showing the first half of the 1994 measure- 
ments of Tau. The curves drawn represent the 13-frequency fit 
found in a later section. 



A journal of the photoelectric measurements used is 
given in Table 1. 

All observatories used the same two comparison stars as 
during the previous DSN campaigns of 9^ Tau, viz., HR 1422 
and HR 1428. No evidence for any variability of these com- 
parison stars was found. We could combine the V and y 
data, because the y magnitudes are defined to be equal to 
the V magnitudes. After the standard photometric reduc- 
tions relative to the two comparison stars, the data from 
the different observatories needed to be combined. For the 
initial multifrequency analyses, this was done by subtracting 
the average relative magnitude of 9^ Tau at each observa- 
tory and then combining the data. Since the nmltifroquency 
analyses utilized the fitting of sine curves to the data, the 
zero-points could then be further adjusted from the values 
of zero-point of these fits. The times of observation were 
converted to Heliocentric .lulian Date (HJD). We also cal- 
culated the light time corrections for the orbital motion of 
the primary stellar component (see Fig. 6 of Paper I). Dur- 
ing this short campaign of 25d, the light time correction is 
very small and practically constant. Therefore, we have not 
applied such a correction to the present data. The extinction 
coefficients were derived separately for each night by using 
the two comparison stars. For a few short data sets, average 
extinction coefficients had to be applied. 
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Table 1. Journal of photoelectric measurements of 6'^ Tau obtained during the 12th Delta Scuti Network campaign 
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Lo: Lowell Observatory, USA; McD: McDonald Observatory, USA; SNO: Sierra Nevada Observatory, Spain 
XL: Xing-Long Observatory, China 



3 METHODS USED TO DETERMINE 

MULTIPLE FREQUENCIES FROM THE 
LIGHT CURVES 

The pulsation frequency analyses wore performed with a 
package of computer programs with single-frequency and 
multiple-frequency techniques (PERIOD98, Sperl 1998). 
These programs utilize Fourier as well as multiple-least- 
squarcs algorithms. The latter technique fits a number of 
simultaneous sinusoidal variations to the observed light vari- 
ability and does not rely on prewhitening. For the purposes 
of presentation and initial searches, however, prewhitening 
is required if the low-amplitude modes are to be seen. There- 
fore, in the presentation of the results (sec below), the var- 
ious power spectra are presented as a series of panels, each 
with additional frequencies removed relative to the panel 
above. 

Two observatories also provided measurements in a sec- 
ond filter, viz., the Stromgren v filter. In this filter, the am- 
plitudes of 0^ Tau are higher by about 30% relative those 
in the y filter. In power spectra, the location of the peaks 
in frequency are not affected by mixing data with different 
amplitudes. Consequently, for the power spectra we have 
used the data from all filters. The slight phase shift up to 
10° in the light curves of 5 Scuti stars between the different 
filters was ignored, since numerical simulations show that 
for 5 Scuti stars the power spectra are not affected by such 
small shifts. 

However, for prewhitening, the data from the two fil- 
ters had to be treated separately. The multipcriodic fitting 
of sinusoids to the data relies on the minimization of the 
residuals between the fit and the observations and requires 
correct amplitudes. This problem was solved by computing 
separate solutions for the two colors, which could then be 
prewhitened. 

One of the most important questions in the examination 
of multiperiodicity concerns the decision as to which of the 
detected peaks in the power spectrum can be regarded as 
variability intrinsic to the star. Due to the presence of non- 
random errors in photometric observations and because of 
observing gaps, the predictions of standard statistical false- 
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Figure 2. Light curves showing the second half of the 1994 mea- 
surements of 6'^ Tau. 



Figures 1 and 2 show the observed light variations of 
9^ Tau together with the 13-frequency fit derived in later 
sections. The new data will be made available in the Com- 
munications in Asteroseismology. 
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Table 2. The photometric frequency spectrum of 6'^ Tau during 
1994 
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Figure 3. Power spectra of the 1994 photometry of 9^ Tau. 
The top panel shows the spectral window, while the other panels 
present the power spectra before and after prewhitening a given 
number of frequencies. See text for a discussion of significance 
levels. 



alarm tests give answers which are considered by us to be 
overly optimistic. In a previous paper (Brcgcr ct al. 1993) 
we have argued that a ratio of amplitude signal/noise = 4.0 
provides a useful criterion for judging the reality of a peaJf. 
This corresponds to a power signal/noise ratio of 12.6. Sub- 
sequent analyses comparing independent data sets have con- 
firmed that this criterion is an excellent predictor of intrinsic 
vs. possible noise peaks, as long as it is not applied to very 
small data sets or at low frequencies, where the errors of 
measurement arc far from random. In the present study, the 
noise was calculated by averaging the amplitudes (oversam- 
plcd by a factor of 20) over 5 cd~^ regions centered around 
the frequency under consideration. 



4 THE MAIN PULSATION FREQUENCIES OF 
TAU DURING 1994 

The spectral window of the 1994 data (Fig. 3) is quite clean 
because the data were collected on three continents, thereby 
avoiding serious regular day-time observing gaps. Neverthe- 
less, the 1 cd"'^ aliasing should still be kept in mind, since 
modes with very small amplitudes might hide in the alias 
peaks. 

The power spectrum of 0^ Tau was computed from fre- 
quency values of zero to the Nyquist frequency. The highest 
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power levels were found in the 10 to 16 cd~^ region. An 
additional band of power was found in the 26 to 30 cd^^ 
region and will be examined in a later section. The power 
spectra in the 8 to 16 cd~^ region are presented in Fig. 3. 
The second panel from the top shows the dominant three 
frequencies, which have been found in all other photometric 
campaigns as well. The next panel shows the power spec- 
trum after prewhitening of the main three frequencies. Three 
additional peaks are detected without any problems. 

The power spectrum of the data after prewhitening six 
frequencies (Data - 6f) presents a small problem: the domi- 
nant peak may be double. The peak at 11.770 cd~^, called 
fr, is strong and its frequency easily determined. Removal 
of this mode leaves another peak at 11.730 cd~^ at a lower 
amplitude. The frequency separation of 0.04 cd^^ is at the 
limit of frequency resolution of the present data set. This 
prevents us from applying the powerful phasing-amplitude 
test to distinguish between two independent close frequen- 
cies and the artifacts caused by a single frequency with vari- 
able amplitudes (see Breger & Bischof 2002) . Fortunately, it 
is not necessary here to prove the existence of a frequency 
doublet: both frequencies have been seen before in the WIRE 
data (Paper III), while the weaker component in the 1994 
data was already detected in the 1986 data (Paper II). We 
therefore accept the reality of the two close frequencies, fy 
and fs. 

Three additional frequencies are found in the 1994 data. 
Further peaks are below the statistical threshold (bottom 
panel of Fig. 3). 

The adopted frequencies and amplitudes are shown in 
Table 2. The uncertainties in amplitude were calculated from 
the standard formulae given in Breger et al. (1999), based 
on the assumption of random observational errors and no 
aliasing. The uncertainties of the amplitudes of fr and fs 
are larger because these two modes are very close in fre- 
quency. Because the separation of two frequencies is at the 
resolution limit of the present data set, the computed size 
of the amplitude of one mode is somewhat affected by the 
presence of the other mode (the so-called pumping effect). 
These amplitudes are therefore marked with a colon. 

The new solution fits the observed data well, as shown 
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in Figs. 1 and 2. The good fit is also demonstrated by the 
value of the average residual per single measurement of ±2.7 
mmag in v and ±3.2 mmag in y. 



5 COMPARISON WITH PREVIOUS STUDIES 

There now exist three large photometric studies of 6^ Tau: 
the present study based on 1994 multisite data, the re- 
analysis of the 1982-1986 multisite campaigns (with empha- 
sis on the 1986 data, Paper II), and the WIRE satellite study 
from space (Paper III). Each photometric study detected 
between 10 and 12 frequencies. The excellent agreement be- 
tween the frequencies found in these independent studies is 
remarkable. 

5.1 The main modes: fi to fr,, fg to fn 

These 9 frequencies were detected in all three studies, but 
with variable amplitudes. While during the years 1982, 1983 

and 1986, the mode at 13.23 cd~^ is dominant in each of 
these years, during 2002 the 13.69 cd^^ frequency became 
dominant. 

5.2 The mode at 10.86 cd^: fe 

The mode at 10.86 cd~^ is present with amplitudes in excess 
of 1.0 mmag in both the 1994 and 2000 data, but appears 
to be absent {y amplitude < 0.3 mmag) during 1986. 

5.3 The 11.73/11.77 cd"^ doublet: fr and fs 

The frequency separation of this close frequency pair is near 
the frequency resolution given by the 25d length of the 1994 
data set as well as that of the 1986 and the WIRE results. In 
the WIRE data, the two frequencies are both present with 
similar amplitudes. During 1994, 11.77 cd^^ is dominant but 
the other peak is also visible. In the 1986 data, we can only 
see the other, lower-frequency mode. The 1982/3 data offer 
higher frequency resolution, but are not useful for the study 
of the pair because of a higher noise level. 

We conclude that both modes arc probably present. The 
strong amplitude variability of these modes may be intrinsic 
to the star or an artefact caused by poor frequency resolu- 
tion. 

5.4 Additional frequencies at 12.70 and 13.81 cd^ 

Two additional frequencies are seen at statistically signifi- 
cant levels in only one of the three data sets: 

(i) The 1986 data reveal an additional mode at 
13.81 cd~^, which has not been seen before. The power spec- 
trum of the WIRE residuals shows a peak at 13.83 cd~^ (see 
Fig. 6 of Paper III) and is therefore probably also present 
in that data set. 

(ii) A mode at 12.70 cd"^ was also found by WIRE, 
which is not affected by 1 cd~^ aliasing and is therefore not 
affected by another mode near 11.73 cd^^. A formal solution 
with the 1994 data including this frequency yields a minis- 
cule amplitude of 0.3 mmag. The mode may be present in the 
1986 data, because in the power spectrum a small peak at 



an amplitude level of ~ 0.5 mmag is seen at that frequency 
value, but below the level of significance. Numerical simula- 
tions suggest that the slight 1 cd^^ aliasing still present in 
the multisite ground-based campaign data is not responsi- 
ble for the computed small (or zero) amplitudes in 1986 and 
1994. 



5.5 Amplitude ratios and phase difFerences 

The variability measurements from Lowell and Sierra 
Nevada Observatories were obtained with the two v and y 

filters of the uvbyP system in order to compare amplitude 
ratios and phase differences in different colors. Especially the 
phase differences for different pulsation modes allow a dis- 
crimination between different I values, but extensive mea- 
surements of extremely high photometric accuracy are re- 
quired. The present study, regrettably, is not large enough 
for reliable mode identifications. In fact, the formal uncer- 
tainties in the phase differences are less than ±3.0° for only 
three modes: fi, f2 and ii. All three modes give similar am- 
plitude ratios with an average value of Av/ Ay = 1.25 ± 0.05. 
For these modes the phase differences, 4>v — <t>y^ s-rc -7° with 
uncertainties of ±2° for the first two modes and ±3° for U,- 



6 THE HIGH-FREQUENCY REGION 

An additional frequency region with power in excess of the 
expected noise occurs in the 26 to 30 cd~^ region. This is 
shown for the 1994 data in the top panel of Fig. 4. Several 
promising peaks exist, but do not reach the amplitudes re- 
quired for a statistically certain discovery of a new pulsation 
frequency, viz., amplitude signal/noise ratios > 4.0. 

The statistics can be improved considerably by exam- 
ining other photometric data in the literature, viz. the data 
obtained during 1982, 1983 and 1986 (see Paper I). For the 
1982/1983 and 1986 photometry, we have prewliitened the 
variability in the low-frequency domain and computed power 
spectra in the high-frequency region. This is also shown in 
Fig. 4. The peak at 26.2 cd~^ is present in all data sets. 
In spite of the small amplitude of 0.6 mmag, the detections 
are statistically significant in two of the three data sets. An 
additional peak at 26.7 cd~^ is also statistically significant 
in the 1982-1986 data, but not seen in 1994. On the other 
hand, the most promising peak in 1994, 28.9 cd~^, is not 
seen in the earlier data. Since its detection is not statisti- 
cally significant in any (including a combined 1982 - 1994) 
data set, it is not adopted by us. 

It is interesting that the 26.2 cd~^ variability was also 
seen independently by the WIRE satellite: Poretti et al. in 
Paper III ascribe this to a spurious term caused by the duty 
cycle and experimental errors. In view of the present de- 
tection of the same frequency in eaxth-based data, we can 
rule out an experimental origin and conclude that these high 
frequencies originate in 9^ Tau. 

Table 3 lists the amplitudes for the two significant 
modes. 
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Figure 4. Power spectra in tlie liigli-frequency region. The fre- 
quency at 26.2 cd~^ is present in all data sets, while 26.7 cd~^ 
is present in both the 1986 and 1982/3 data sets. 



Table 3. Variability of 0'^ Tau at high frequencies 
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Amplitudes^ 
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1986 


1994 V 


1994 y 




cd-i 




mmag 


mmag 


mmag 


mmag 








±0.1 


±0.1 


±0.1 


±0.1 


fl2 


26.184 


303.1 


0.6 


0.6 


0.5 


0.4 


fl3 


26.732 


309.4 


0.5 


0.5 


0.1 


0.1 



^ Note that an amplitude of 0.1 mmag does not represent a sta- 
tistically significant detection. 



6.1 Could the high frequencies be combination or 
harmonic frequencies? 

S Scuti stars often show combination frequencies, fi±fj, as 
well as slightly nonsinusoidal light curves leading to 2fi peaks 
in the power spectrum. These peaks are, of course, situ- 
ated at high frequencies and axe characterized by amplitudes 
smaller than a factor of ten (or more) relative to those of 
fi. Since 9^ Tau has its main variability in the 11-15 cd~^ 
range, combination frequencies/2f terms would be found in 
the 22-30 cd~^ region and therefore provide a potential ex- 
planation. 

The frequencies of the potential high-frequency candi- 
dates can be easily calculated from fi±fj and 2fi. Such have 
been detected even in low-amplitude 6 Scuti stars such as 
4 CVn (Breger et al. 1999) and XX Pyx (Handler et al. 1996, 
2000). 

The value of the main high frequency, fi2 at 
26.184 cd~^, cannot be identified with the frequency of a 
combination of any two modes listed in Table 1. We conclude 
that it is not a combination frequency. On the other hand, 
fi3 at 26.732 cd~^, is in the vicinity of two potential combi- 



nations, fi+f3 = 26.714 and f4+f5 = 26.713 cd"\ Although 
the 1982-1986 data have an excellent frequency resolution 
due to the 1000+ day span, the data contain long gaps, 
leading to the possibility of spectral leakage. We can see 
smaller peaks near 26.714 cd"\ so that some of the power 
at 26.732 cd~^ may be spectral leakage from combination 
frequencies. The important question concerns the physical 
origin of the combination terms, e.g., are these combination 
frequencies just combination terms or are new modes excited 
by resonance at frequencies near those of the combination 
frequencies? Garrido & Rodriguez (1996), in their study of 
combination frequencies, point out that the phasing of these 
combination terms might allow us to distinguish between the 
two possibilities. However, for low-amplitude 6 Scuti stars, 
accurate phasing of combination terms are observationally 
difficult to determine and the question cannot be examined 
further at this stage. 

Observations of S Scuti stars (and other pulsating vari- 
ables) with known combination frequencies show that the 
amplitudes of the combination terms are a function of the 
amplitudes of the two frequencies involved in the combi- 
nation terms. We can estimate the size of the combination 
terms, fi+fj and 2fi from the Garrido & Rodriguez analysis 
as well as the observations of XX Pyx, 4 CVn and BI CMi 
(Breger et al. 2002) and calculate expected amplitudes small 
than 0.1 mmag. This is an order of magnitude smaller than 
observed. The amplitude argument, therefore, does not sup- 
port an identification of the two frequencies as combination 
modes. 

We conclude that the main peak, fi2 cannot be identi- 
fied as a combination frequency, while such an explanation 
is unlikely for fis. 

6.2 Pulsation in a 5 Scuti binary companion 

6^ Tau is a 140.728d binary system with two components 
of similar temperature. Peterson, Stefanik & Latham (1993) 
have shown that the secondary is fainter by 1.10 mag and 
still on the main sequence. Both stars are situated inside the 
instability strip. It was shown in Paper I that the dominant 
pulsation modes originate in the evolved primary component 
because the predicted orbital light-time effects for the pri- 
mary match the observed shifts in the light curves of up to 
several minutes. This confirms an expectation from the val- 
ues of the frequencies (10 - 15 cd~^) which are incompatible 
with those of main-sequence 5 Scuti stars. 

What would be the expected frequencies of pulsation of 
the main-sequence secondary? Although S Scuti stars pul- 
sate with many simultaneously excited modes of pulsation, 
the average time-scale obeys a period-luminosity relation of 
the form, M„ ~ -3.05 log P (see Breger 2000). We can now 
apply the known magnitude difference between the compo- 
nents to provide a rough estimate of the frequencies expected 
for the secondary. The magnitude difference of 1.10 mag 
shifts the 11-15 cd~^ range of the primary to 25-33 cd~^, 
which agrees with the band of power seen in the power spec- 
trum. 

The light of the secondary is heavily diluted by the light 
from the primary, which is 1.1 mag brighter than the sec- 
ondary. Consequently, the amplitudes of pulsation of the sec- 
ondary are reduced to one fourth of their intrinsic values. If 
our interpretation is correct, then the secondary would have 
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V amplitudes of about 0.5 mmag in combined light, or 2 
mmag in the star alone. The hypothesis that these modes 
come from the secondary can be tested: the orbital motions 
produce light-time shifts of several minutes. Regrettably, the 
observed amplitudes are too small to obtain the required ac- 
curacy in the phasing at different orbital times. 



7 PULSATION MODELS FOR 9^ TAU 

In order to examine the nature of the pulsations of 0^ Tau 
in more detail, the physical parameters for the two stars 
forming the binary system need to be determined. A recent 
summary of different determinations of global parameters of 
0^ Tau can be found in de Bruijne et al. (2001). Torres et 
al. (1997) used the orbital parallax to derive absolute visual 
magnitudes, Mv, of 0.37 and 1.47 for the two components. 
Similar values, 0.48 and 1.58, are obtained by de Bruijne et 
al. (2001) and Lebreton at al. (2001) from Hipparcos dynam- 
ical parallax. Both stars have similar temperatures and its 
value can be obtained from the Moon & Dworetsky (1985) 
calibration of uvbyf3 photometry, viz., Teff = 7950 K. This 
value is lower than the value of 8200 K adopted by Breger et 
al. (1987) as well as others, since the older model-atmosphere 
calibrations of uvby/3 tended to overestimate the tempera- 
tures of A stars. Using a calibration of the UBV photometry, 
de Bruijne et al. (2001) derive Tcff = 7980 K for the primary 
and Teff — 8230 K for the secondary component 

The Moon & Dworetsky (1985) calibration also allows 
us to derive a log g = 3.69 from the uvbyf3 photometry of 
the combined light of the two stars. Correction of the mea- 
sured Balmer discontinuity, ci , for the less evolved secondary 
places our estimate for the primary near log g = 3.55. 

This implies a considerably higher evolutionary status 
than that deduced by Torres et al., who adopted log g = 4.0. 
However, our value is consistent with the position of the star 
in the Hertzsprung-Russell Diagram for the Hyades cluster 
as well as the known luminosity (see above). We can also 
derive the log g value from the known mass, luminosity and 
temperature: for a mass of 2.42 solar masses for the primary 
(Torres et al.) a log g value of 3.63 is derived, in agreement 
with the value found above from uvbyP photometry. 

Finally, this as well as similar studies neglect the effect 
of rotation and aspect, which affect the choice of the appro- 
priate values for temperature, luminosity and gravity (e.g., 
see Perez Hernandez et al. 1999). 



7.1 Method of computation 

The method of computation of stellar evolution and pulsa- 
tion was the same as in our other studies of 5 Scuti stars (see 
Pamyatnykh 2000 and references therein). In particular, we 
used the latest version of the OPAL opacities (Iglesias & 
Rogers 1996) supplemented with the low-temperature data 
of Alexander & Ferguson (1994). Also, we used an updated 
version of the OPAL equation of state (Rogers et al. 1996), 
viz., version EOS2001, which was copied from the OPAL 
data library^. 



I |ht.t.p: //wM-phva . 11 nl . gov/RRSRarch/OPAL/opal . html 

f tp: //www-phys . Ilnl.gov/pub/opal/eos2001/ 



The computations were performed starting with chemi- 
cally uniform models on the ZAMS, assuming an initial hy- 
drogen abundance X = 0.716, a helium abundance Y = 0.26 
and a heavy element abundance Z — 0.024. These values 
correspond to recent estimates of the chemical composition 
of the Hyades (Ferryman et al. 1998, Lebreton et al. 2001). 
We also tested models of somewhat different chemical com- 
position, with slightly higher metallicity or with higher he- 
lium content {X = 0.723, Y = 0.25, Z = 0.027; X = 0.686, 
Y = 0.29, Z = 0.024). 

We computed models with and without overshooting 
from the convective core. In the first case, the overshooting 
distance, dover, was chosen to be 0.2 Hp, where Hp is the lo- 
cal pressure scale height at the edge of the convective core. 
Examples of evolutionary tracks for S Scuti models com- 
puted with and without overshooting are given by Breger & 
Pamyatnykh (1998) and Pamyatnykh (2000). In the stellar 
envelope the standard mixing-length theory of convection 
with the mixing-length parameter a = 1.6 was used. This 
value was also used by Ferryman et al. (1998) and Lebre- 
ton et al. (2001), who derived a from the calibration of the 
solar-model radius. 

Uniform (solid-body) stellar rotation and conservation 
of global angular momentum during evolution from the 
ZAMS were assumed for all computations. These assump- 
tions were chosen due to their simplicity. The influence of 
rotation on the evolutionary tracks of S Scuti models was 
demonstrated by Breger & Pamyatnykh (1998) and Pamy- 
atnykh (2000). In most cases the initial equatorial rotational 
velocity on the ZAMS was chosen to be 100 km/s. 

The linear nonadiabatic analysis of low-degree oscilla- 
tions {£ < 2) was performed using the code developed by 
Dziembowski (1977). The effects of slow rotation on oscil- 
lation frequencies were treated up to third order in the ro- 
tational velocity (Dziembowski & Goode 1992, Soufi et al. 
1998). 

The main aim of the present theoretical study was to 
check the instability of models of both stellar components 
of 6^ Tau for the whole range of allowed values of mass or 
gravity, effective temperature and luminosity and to outline 
possible asteroseismological constraints to the models. We 
did not attempt to construct any model which fits quanti- 
tatively the observed and the theoretical frequency spectra 
of excited multiperiodic oscillations. 

As a first step, we tried to fit the observed and theo- 
retical frequency ranges of the unstable modes. The main 
uncertainty of such a study is the unsatisfactory description 
of convection in the stellar envelope and its interaction with 
pulsations. The fact that convection may influence the insta- 
bility of stars not only near the Red Edge but almost in the 
whole classical instability strip was discussed by Stellingwerf 
(1984) in connection with the problem of the determination 
of helium abundance in globular clusters from the tempera- 
ture at the Blue Edge of the RR Lyrae instability domain. 
5 Scuti stars are generally hotter than RR Lyrae stars, but 
even in this part of the instability strip convection influences 
the onset of the instability for low overtones of radial and 
nonradial oscillations. Fig. 9 of Pamyatnykh (2000) demon- 
strates that convection signiflcantly displaces the radial fun- 
damental Blue Edge to hotter temperatures due to increased 
driving in the hydrogen ionization zone which is more ex- 
tended than in the radiative models. A part of this contribu- 
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tion to the driving may also be caused by the generally-used 
assumption that the convective flux does not vary during an 
oscillation cycle. This simple assumption is probably incor- 
rect in the hydrogen convection zone. Therefore, the real- 
ity of the additional convective dcstabilization in the outer 
hydrogen zone must be examined by using a nonlocal time- 
dependent treatment of convection. First promising results 
in this direction (see Michel et al. 1999, Houdek 2000) ex- 
ist. Note that the total driving in a 5 Scuti star (with main 
contribution due to the k mechanism operating in the sec- 
ond helium ionization zone) only slightly exceeds the total 
damping. Therefore even small contributions to the driving 
are important. 

Due to these uncertainties wc consider our results as 
preliminary and our models as test models. Nevertheless, 
even with such a simple description of the convection we 
may have a possibility to constrain the stellar parameters 
from instability studies. 

7.2 Results for test models 

In Fig. 5 we show evolutionary tracks for selected models, 
computed by taking overshooting from the stellar convec- 
tive core into account. The models with effective tempera- 
ture Teff = 8000, 7900 and 7800 K, in which we study the 
oscillations, are marked by filled circles. The luminosities of 
the higher and the lower-mass models approximately corre- 
spond to the luminosity of the primary and secondary stellar 
components, respectively. 

For the primary component, we have chosen models at 
the beginning of the post-main-scqucnco (hereafter called 
post-MS) expansion stage, which burn hydrogen in a thick 
shell just after full hydrogen exhaustion in the stellar core. 
This evolutionary stage of the primary was proposed earlier 
by Krolikowska (1992) for a significantly hotter model of 
about 8200 K without overshooting. Moreover, wc studied 
slightly more massive models of similar effective temperature 
and luminosity in the MS stage as well as post-MS models 
without overshooting. 

The secondary component, which is fainter by approxi- 
mately 1.1 mag (which corresponds to the diff'crcncc by 0.44 
in logL between the components) is a main-sequence star 
with hydrogen burning in the convective core. 

In Fig. 5 we also plot the Blue Edges of the instability 
strip taken from Pamyatnykh (2000), which were computed 
for nonrotating models without overshooting for X = 0.70, 
Y = 0.28, Z = 0.02 and with an assumed mixing-length 
parameter a = 1.0. The general Blue Edge is defined as 
the hotter envelope of unstable overtones, from radial mode 
ps (seventh overtone) near the ZAMS to mode (fourth 
overtone) for highest luminosities or smallest gravities in 
the figure (see Fig. 3 in Pamyatnykh 2000). The Blue Edge 
for the radial fundamental mode lies approximately in the 
center of the 5 Scuti instability strip. 

An additional study of the instability along the 2.5 Mq 
evolutionary track shows that the best theoretical general 
Blue Edge for X = 0.716, Y = 0.26, Z = 0.024 will be lo- 
cated very close to the blue edge shown in Fig. 5, because 
the differences in the rotational velocity and in the over- 
shooting efficiency do not infiuence the position of the Blue 
Edges and because the differences in the helium abundance 
are small. Moreover, convection has only a minor influence 
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Figure 5. Evolutionary tracks for indicated values of M/Mq. 
Test models with effective temperature T^g = 8000, 7900 and 
7800 K for the primary and the secondary component of 6'^ Tau 

are marked by filled circles. All models were computed for an 
initial hydrogen abundance X = 0.716, a helium abundance 
Y = 0.26 and a heavy element abundance Z = 0.024. A mixing- 
length parameter of convection a = 1.6 was used. Overshooting 
from the stellar convective core was taken into account with an 
overshooting distance dover = 0.2 ifp. The equatorial rotational 
velocity on the ZAMS was chosen to be 100 km/s. The almost ver- 
tical lines show the general Blue Edge of the instability strip and 
the Blue Edge of the radial fundamental mode, pi, as computed 
for a standard chemical composition (Pamyatnykh 2000). 

on the position of this hot general Blue Edge (see Fig. 9 
in Pamyatnykh 2000) . For the fundamental radial mode the 

best Blue Edge will bo hotter by 0.008 - 0.009 in logTcff. 
This is mainly due to a higher value of the mixing-length 
parameter. 

From Fig. 5 we immediately obtain a strong constraint 
on the possible effective temperature of the primary of 
e"^ Tau. All models with log(L/L0) > 1.67 and log Toff > 
3.907 (Teff > 8070 K) are stable in all modes. A MS model 
of 2.5 Mq on the Blue Edge (Tch = 8070 K) is marginally 
unstable in radial mode pe with the frequency 18.74 cd"'^, 
which is well outside the observed frequency range. We can 
conclude that only significantly cooler models can pulsate 
with the observed frequencies in the 10.8 to 14.6 cd~^ range. 

This conclusion is confirmed by computation of oscilla- 
tions of the selected test models for the primary of Tau. 
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Figure 6. Normalized growtli rates, rj, plotted against frequency, 
/, in test models of the primary comijoncnt of 9^ Tau. Positive 
values of r] correspond to unstable modes. The values of M (in 
solar units) and Tgg are given in each panel. The position of 
the models in the HR diagram is shown in Fig. 5. Models were 
calculated with initial abundances X = 0.716, Y = 0.26, Z = 
0.024 and VrotCZAMS) = 100 km/s. The symbol ps is plotted 
near the corresponding radial overtone. The thick horizontal line 
shows the range of frequencies (10.865 to 14.615 cd~^) observed 
for the primary component of 6^ Tau. 



In Fig. 6, the normalized growth rates of radial and nonra- 
dial modes arc plotted against frequency for all nine higher- 
mass models which are marked in Fig. 5. Only axisymmetric 
modes (m = 0) are shown. The independence of the growth 
rate on the spherical harmonic degree, £, is a typical fea- 
ture of modes excited by the k mechanism. The rotational 
velocities of the models are 81 to 86 km/s. The rotational 
splitting of the modes can extend the frequency range by ap- 
proximately 0.5 c/d on both sides. We can see that excited 
frequencies of the 2.45 Mq model with Tes = 7800 K are 
in excellent agreement with the observed frequency range. 
The frequency range of unstable modes spans three radial 
orders from p4 to pe for radial modes (mode p4 is marginally 
unstable) . 

As was noted already, the results are sensitive to the 
treatment of convection. For example, if we use a mixing- 
length parameter a = 2.0 instead of a = 1.6, the fre- 
quency range of unstable modes for 2.45 Mq model with 
Toff — 7800 K is extended by 1 c/d on both sides. Moreover, 
our assumption about the unperturbed convcctive flux dur- 
ing an oscillation cycle is not fulfilled inside the hydrogen 
convective zone and may result in artificial additional driv- 
ing in this zone. Therefore, these preliminary results must 
be considered with caution. 

Similar results were obtained for the models without 
overshooting, the best fitting is achieved in this case for 
2.50 Mq model with T^a = 7800 K. Also, for slightly more 
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Figure 7. The same as in Fig. 6 but for test models of the 
secondary component. Short vertical lines mark the observed fre- 
quencies 26.184 and 26.732 cA'^ . 



massive MS models with overshooting, we obtained a good 
agreement between the observed and the theoretical fre- 
quency ranges. The parameters of some models are given 
in Table 4 below. 

In Fig. 7 we show the normalized growth rates in test 
models of the secondary component. These models are also 
marked in Fig. 5. As for the primary, we used models with 
effective temperatures from 7800 to 8000 K. The frequency 
range of unstable modes is much wider for these less lumi- 
nous models: it spans 5 to 7 radial orders from p2 to ps- 
The observed frequencies lie near or between radial modes 
P5 to pe, which are unstable in the whole range in effective 
temperature. 

In Table 4 we summarize parameters of the best test 
models. We can see that both MS and post-MS models for 
the primary with T^s = 7800 K can fit the observed fre- 
quency range. The values of the surface gravity and lumi- 
nosity of all three models are similar. In the computation 
of the evolution of the MS model PI, the initial equatorial 
rotational velocity on the ZAMS was chosen to be 80 km/s. 
(Rotation does not affect our theoretical results significantly, 
but it may influence the photometric calibrations.) The to- 
tal number of unstable modes of £ < 2 is equal to 37 in the 
MS model PI, 77 and 56 in the post-MS models P2 and P3, 
respectively. 

Note that the post-MS evolution of models for the pri- 
mary is much faster than the MS evolution of the less mas- 
sive secondary. For example, in Fig. 5 all models of 2.45 Mq 
with Toff = 8000 K, 7900 K and 7800 K (the last model is just 
model P3) have the same age 695 Myr. These models corre- 
spond to stars which are slightly older than the Hyades (650 
Myr, Lebreton et al. 2001). On the other hand, the models 
without overshooting may correspond to stars younger than 
the Hyades: the age of 2.50 Mq models is 595 Myr. The best 
model for the primary may therefore be one with slightly 
weaJjer overshooting than that of model P3. 

The inferred age of the secondary star can be adjusted 
to the age of the primary as well as to that of the Hyades 
cluster. This is possible because of the slower main-sequence 
evolution of the secondary and due to the fact that corre- 
sponding models of different effective temperatures flt the 
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Table 4. Parameters of the best test models of 0-^ Tau which fit the observed frequency ranges. Models of 2.65, 2.50 and 2.45 Mq 
correspond to the primary component, models of 2.10 Mq correspond to the secondary component. The symbols have their usual 
meanings (see text). The last column gives the frequency range, in cycles per day, of unstable modes for each model. 



Model 


M/Mq 


Ev. Stage 


Age [Myr] 


Teff 


logL 


log g 


Vrot 




Frequency range 


PI 


2.65 


MS 


535 


7800 


1.792 


3.58 


60 


0.2 


10.7- 15.4 


P2 


2.50 


post-MS 


595 


7800 


1.780 


3.56 


83 


0.0 


10.5 - 14.6 


P3 


2.45 


post-MS 


695 


7800 


1.762 


3.57 


84 


0.2 


10.6 - 14.9 


SI 


2.10 


MS 


667 


7900 


1.319 


3.97 


91 


0.2 


13.0 - 37 


S2 


2.10 


MS 


704 


7800 


1.325 


3.94 


90 


0.2 


16.5 - 38 



observed frequencies quite well. As an example, models of 
2.10 Mq with Tefi = 8000 K, 7900 K and 7800 K have the 
ages 625, 667 and 704 Myr. 

Moreover, we tested some models of different chem- 
ical composition. If we choose the estimated values for 
the helium content (Lebreton et al. 2001) and metallicity 
(Lastennet et al. 1999), the chemical compositions become 
X = 0.723, Y = 0.25, and Z = 0.027. The results for 
these models arc very similar to those computed above for 
X = 0.716, Y = 0.26, Z = 0.024. 

Models with X = 0.686, Y = 0.29, Z = 0.024 are some- 
what more luminous for a given mass, e. g., the track of 
2.45 Mq almost coincides with the track of 2.55 Mq com- 
puted with X = 0.716, Y = 0.26, Z = 0.024. The models 
are also pulsationally more unstable and the range of un- 
stable modes is somewhat wider due to higher helium abun- 
dance. As a consequence, the MS model of 2.45 Mq with 
TefF = 7900 K can be unstable in the observed frequency 
range and also at slightly higher frequencies. The luminosity 
of this model, log(L/Lo) = 1.71, agrees with the estimates 
from the orbital and dynamical parallaxes. 



8 CONCLUSION 

The DSN campaign 12, carried out during 1994 at four ob- 
servatories, resulted in 152 hours of high-precision photom- 
etry. 13 frequencies of pulsation were derived from the data. 
These frequencies confirm the results from previous Earth- 
based (1982-1986) as well as satellite (2000) photometry, al- 
though amplitude variability on a time scale of several years 
is present. 

While the detected modes in the 10-15 cd~^ range prob- 
ably all originate in the primary, the higher frequencies with 
small amplitudes probably originate in the 5 Scuti-type of 
variability of the main-sequence companion. 

We constructed evolutionary models both for the pri- 
mary and the secondary components and tested them for 
the instability against radial and nonradial oscillations. 

The best fit of the theoretical and observed frequency 
ranges is achieved for models with T^s ~ 7800 K or slightly 
higher, in agreement with photometric calibrations. The in- 
stability spans two or three radial orders from p4 to pe for 
radial modes. Post-MS models with or without overshooting 
are preferable for the primary, but MS models with over- 
shooting are also possible. 

For the less luminous secondary component the insta- 
bility range is wider and spans 5 to 7 radial orders from p2 



to ps- Models of different effective temperature are unstable 
at observed frequencies 26.184 and 26.732 cd~^ which lie 
around radial modes ps to pe- 

The main uncertainties of the results arc caused by our 
crude treatment of the convcctivc flux. We neglected the 
Lagrangian variation of the convcctive flux during an oscil- 
lation cycle. This assumption is not fulfilled in the hydrogen 
convective zone, which may lead to an incorrect estimate 
of the small contribution of this zone to the driving. To 
prove the results quantitatively, it is necessary to use a reli- 
able theory of non-local time-dependent convection. We note 
also the nonlocal results by Kupka & Montgomery (2002), 
who point out the necessity to use very different values of 
the mixing-length parameter in the hydrogen and helium 
convection zones, if we still use local mixing-length theory. 
According to nonlocal studies of A-star envelopes, it is nec- 
essary to use a small value of this parameter in the hydrogen 
zone and a significantly higher value in the deeper ficlium 
zone. We plan to perform corresponding tests for the 9^ Tau 
models. However, even with the present calculations and re- 
sults we probably can exclude modes involving the second 
radial overtone (mode ps) and lower-order modes, as well 
as all hot models of the primary with T^g > 8000 K. The 
reason is that all these models are stable in the observed 
frequency range of 10.8 to 14.6 cd~^. 
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